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Abstract 

We obtain the approximate solutions for the steady temperature profiles of materials with a 
temperature-dependent thermal absorptivity inside a microannulus with wavy-rough surfaces con- 
sidering a quasilinear partial differential equation by the boundary perturbation approach. Our nu- 
merical results show that the critical Frank-Kamanestkii parameter depends on the small-amplitude 
wavy-roughness for specific power-law of thermal absorptivity. 
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I. INTRODUCTION 



Orbiting spacecraft are subjected to forces such as drag and solar pressure that have to 
be compensated for in order to maintain the desired orbit, and they may also have to be 
transferred from the initial orbit where they are placed by the launcher to their final orbit. 
The orientations and orbits of spacecraft are controlled by low-thrust propulsion systems. 
The propulsion control systems can be broadly classified as suitable for coarse and fine con- 
trol. Hitherto, the major effort has been expended on the development of coarse control 
systems. In the 1960s, simple cold gas propulsion systems were employed and over the 
years more efficient systems were developed such as monopropellant hydrazine fuel, storable 
bipropellant mixtures providing more efficient combustion and, quite recently, new type of 
propulsion systems to provide higher propulsive impulse and higher performance to cope 
with increased spacecraft masses and longer satellite lifetimes [1-2]. The crucial issue which 
is related to the important components of many space systems is the thruster system that 
shall generate thrust in a well-defined direction. One means for thrust generation is solid- 
state micro-thrusters which requires the safety storage of fuels [3]. 

Meanwhile, thermal instability such as the catastrophic phenomenon of thermal runaway in 
which a slight change of external heating (e.g., microwave) power causes the temperature 
to increase rapidly to the melting point of the material (e.g., ceramic) have been studied 
for simple geometry cases [4-10]. Up to now, the basic understanding of the (microwave) 
heating process still remains somewhat empirical and speculative due to its highly nonlinear 
character. The mathematical description of this process is accordingly fraught with non- 
linearities (not to mention the external heating) : The heat equation is highly temperature 
dependent at the temperatures required for sintering. Moreover, the thermal boundary con- 
ditions must take into account both convective and radiation heat loss. The result is a highly 
nonlinear initial-boundary value problem. Even for simple geometry problems, because of 
the highly nonlinearity characteristic, numerical approaches were often adopted to solve the 
corresponding equations and boundary conditions [4,9-10]. 

Relevant approaches have also a wide field of applications in combustion theory or for safety 
of storage of materials capable of exothermic chemical reaction [11]. In this short paper, 
we shall consider the above mentioned problems with simple geometry (say, an annulus) 
in microdomain. However, real surfaces are rough at the micro- or even at the meso-scale 
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[12-14]. We presume the roughness to be wavy-like in transverse direction along the outer 
and inner surfaces of the annulus which thus make the corresponding boundary conditions 
highly nonlinear. With this, we adopt the boundary perturbation approach [12-14] to handle 
the mathematically complicated boundary value problems. Our numerical results show that 
for certain power-law of thermal absorptivity the small amplitude wavy-roughness could 
influence the steady thermal behavior inside a microannulus significantly. 

II. MATHEMATICAL FORMULATIONS 

Considering the forced heat equation 



where T is the temperature, v is the thermal conductivity, 7 is the thermal absorptivity, 
p is the density that is usually assumed constant, c p is the specific heat and is the heat- 
ing source. Here, the thermal absorptivity depends on the square of the external source 
(amplitude). Assuming constant specific heat and that the heating source (amplitude) is 
temperature independent, the appropriate form of above equation becomes 



where might depend on the spatial coordinates. 
A. Smooth Cylinder Geometry 

Without loss of the generality, we consider a class of solutions for the nonlinear reaction- 
diffusion equations 



with the geometry of an infinite cylinder (x = r with r being the radius) and m = —1. Here, 
S is the Frank-Kamanestkii parameter [15] which incorporates the chemical properties of the 
combustible material, the temperature of assembly as well as its geometrical dimensions. 
K(x) = x^ and fj, < 0, fj, is a constant. The critical value 5 cr defines criticality and provides, 
for example, a criterion for safety of storage of materials capable of exothermic chemical 
reaction [11,15]. Thus, 5 cr is important in characterizing the thermal stability properties of 
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the material under consideration. 

For the smooth geometry case, we consider the steady state thermal explosions between two 
infinite cylindrical material undergoing exothermic reaction by using the change of variable 
given by x = exp(y), T n = exp(C/) and = U + (2 + fx)y and we can obtain that closed-form 
solutions which are available for all real values of fx [15]. In fact, equation (2) becomes 

d 2 Q 



dy 2 



+ Sne e = (3) 



and the solutions are 



e & = v or T» = CoX ~ 2 ~" v (4) 

cosh 2 (^f^y-D ) cosh 2 (^f^\nx-D ) 

where Co, Dq are integration constants. 

Note that for the pure cylindrical case (x = r = 0), once fx > —1, we must make sure the 
solution to be well defined considering the limit of x — > : 

T n = ^£ . (5) 



2n5x 2 +v cosh 2 In x - D ^ 



with 

_ (2 + /i) 2 

The boundary conditions for a smooth annulus could be 

T n (l) = 1 and T n (s) = 1, (6) 

where s = ri/r 2 (cf. Fig. 1) and x = r/r 2 . These could help us to fix Cq and D . Under 
this situation, we have 

J^—lns cosh(- A,) = Do + cosh -1 [s~^ cosh(-A))] , (7) 

V ~ 

where fx < — 1. 



B. Wavy- Rough Annulus Geometry 

We start to consider a steady state heat flow in a wavy-rough microannulus of r 2 
(in mean-averaged outer radius) with the outer wall being a fixed wavy-rough surface : 
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r = r 2 + esm(kO) and T\ (in mean-averaged inner radius) with the inner wall being a fixed 
wavy-rough surface : r = ri+esin(k9+(3), where e is the amplitude of the (wavy) roughness, 
and the wave number : k = 2ir/L (L is the wave length), (3 is the phase shift. The schematic 
is illustrated in Fig. 1. 

With the boundary perturbation approach from Refs. [12-14], we shall derive the temper- 
ature field along the wavy-rough microannuli considering dimensionless values. We firstly 
select L a = r 2 to be the characteristic length scale and set 



After this, for simplicity, we drop all the primes. It means, now, r, R Q , Ri, and e become 
dimensionless. The walls are prescribed as r = R Q + esin(/c6>), r = Ri + esin(/c# + (5) and 
the presumed steady state heat flow is along the ^-direction (microannulus-axis direction). 
Let T be expanded in e : 



where the subscript means the partial differentiation (say, T r = dT/dr). 

The corresponding equation, for simplicity (as the wavy roughness being imposed) as we 

select firstly n — m + 1, becomes 



r = 



r/L a , R = r 2 /L a , R i = r 1 /L a , e' = e/L a . 



T = T + eTi + e 2 T 2 + • • • , 



(8) 



and on the boundary, we expand T(r + edr, 6{= 9 )) into 



T(r, 9)\ {ro+edrM = T(r , 9) + e[drT r (r , 9)} + e 2 [—T rr (r , 9)} + ■ ■ ■ 




(9) 



or 




(10) 



with = T m+1 . 



Now, we can set rj = r^ 2+fM ^ 2 and then obtain 

<9 2 d& A{m + 1)5 



fe = o, 



(11) 



drj 2 r/drj (2 + /i) 2 



Above equation has the zeroth order solution in terms of Bessel functions 



Q = A J (Kr ] )+B Y Q (Krj), 



(12) 
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once we consider 6 = 6 + e©i + ■ • • [16]. Here, A , B are integration constants. K = 
2<t/(2 + n) and a = [(to + 1)5] 1/2 . Thus, we obtain, for /j ^ -2, 

G = n = To" +I = + W rMl*) + Bo y (*!±l)^! r (-,)/^ (13) 

V 2 / V 2 -I - /x / 

and for /x = — 2, 

o = T o n = T o m+1 = A*sin([(m + l)5] 1/2 lnr) + 5* cos(J(to + l)^ 1 / 2 lnr). (14) 

Here, A ,B or A*,£?* could be uniquely determined after imposing the boundary conditions. 
This leads to T = 0o /(m+1) . 

As for the first order correction or solution, with the wavy- roughness effect, we can only 
have (considering T n = (T + eT x + ■ ■ -) n = T n + enT^' 1 + ■■■ with T known), for to = 0, 

+ <El + 1^ + Sr^Ti = (15) 

due to the complicated orientational characteristics (^-dependent). Using the boundary 
conditions from the equation (6), we have 

To| r=s = 1, 2o| r= i = 1, (16) 

^\ r=1 sm(k6)+T 1 \ r=1 = 0, (17) 
or 

and 

s sm(k9 + p)+T 1 \ r=s = 0, (18) 



dr 

Equation (15) can be transformed to 
d 2 R dR , 2 



and solved after setting 



+ + {^—) 2 {5r» +2 - K 2 )R =0, 7\ = te tKe R{0, (19) 



e = r^ 2 , "=^ lf \ K =^~ K > ^- 2 > ( 20 ) 



where 



j r(0 = a 1 j^K) + 5 1 j_ 1? K)- (2i) 



We thus have 



5+7r x V2 + /i / x V2 + /i / 



(22) 
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with n ^ — 2 or 



Ti = ie ! 



Ci cos ({5 - K 2 ) l/2 lnr) + A sin ((5 - K 2 ) 1/2 \nr 



fi = -2. (23) 

Here, Ai,Bi,Ci,Di could be fixed by boundary conditions or equations (16-18). In fact, 
from the boundary condition (17) or the wavy roughness along the annulus, we can identify 
K — k. 

With the first order perturbed solution due to wavy roughness : Ti, we can thus obtain the 
approximate solution for especially m = 



T = T + eTi + • • • , 



where 



or 



T = A* sin (V /2 In r) + B* cos (V /2 In r) , /i = -2, 

with Ti given by equations (22) and (23). To be specific, we have (s = Ri) 

Y Q {R s s^)l 2 ) - Y (R S ) 

Jo(R s ) Y (R s s(^yz) - Y (R s )J (R sS (^y*y 

MRs) - MRssP + M 2 ) 

Jo(R s ) Y (R 8 se+M*) - Y (R s )J (R s s( 2 +^ 2 ) ' 

where R s = 2y / 5/(2 + //), and 

A* = 1 - C ° S ^ lns) , B* = 1. 

sin(v5 Ins) 

Meanwhile, with equations (17) and (18), we can obtain 

~ F t / 2\fl (9 ^„s/ 9 \ „ t / 2v^ 



A 1 = 



— V2 + /i / — V2 + /iyT 



where 



5+7^ K V2 + /i/ 
Fx = -v^ 



V2 + /i 



s (2+M)/2 \ 



/det, 



, 2v^ x , 2^ 



v 2 + ^' 



d = -e^v 7 ^ 2 



2a/5s( 2+ ^/ 2 

a Ti(^ ) + s y 1 



2 + fl 

2v ^ S (2+M)/2 



(24) 

(25) 
(26) 



(27) 
(28) 

(29) 
(30) 



det = J 2 JP^)J-, K f^ls^ 2 ) - J_„ (M); 2 x fM s ( 2+ ,)/ 2 ) 
~ K V2 + / u/ — V2 + /i / —\2 + fiJ — K \2 + fx J 
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III. NUMERICAL RESULTS AND DISCUSSIONS 

We plot T (r) for different fx (= -4, -f , 1) with fixed 5(— 2) in Fig. f . The boundary 
conditions : T (r = s = 0.4) = T (r = f) = f are correctly matched at both walls. Con- 
sidering the contributions from the first-order perturbation Ti(r), we plot the approximate 
temperature profiles : T (r) + eTi(r) with the same —4, — f , f ) in Fig. 2. Here, e = O.f , 
P = 7r/4 and fc = 8. We can observe that, for \i = —4, the effect of first-order perturbation 
is quite large compared to other cases of /i. 

To consider the effects of roughness, we could observe that, as evidenced from the first-order 
correction or approximation : T ± , the Frank-Kamanestkii parameter (5) depends on the wave 
number of wavy- roughness (cf. Eq. (23)). Thus, the critical Frank-Kamanestkii 5 cr depends 
on the square of the wave number k (or k cr for certain critical ones). However, above results 
are only for cases of \i = — 2 (or K{x) = x~ 2 which is defined in Eq. (2))! 

IV. CONCLUSION 

We adopt the boundary perturbation approach to obtain the approximate solutions for 
the steady temperature profiles of materials with a temperature-dependent thermal absorp- 
tivity inside a microannulus with wavy-rough surfaces considering a quasilinear partial dif- 
ferential equation. Our results show that the critical Frank-Kamanestkii parameter strongly 
depends on the wavy-roughness for specific power-law of thermal absorptivity. Acknowledge- 
ment. The first author stayed at the Chern Shiing-Shen Institute of Mathematics, Nankai 
University for one month. Thus the first author should thank their hospitality for their 
Visiting-Scholar Program. 
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Fig. 1. Approximate zeroth order temperature profiles : ?b(r) for different fi (= —4,-1,1) 
with 5 = 2. The boundary conditions are satisfied at both sides (Tq(s = 0.4) = Tq(1) = 1). 



10 



2.5r- 
2 

1.5- 



9. ? ■* 9- 



O T -|l= -4 

V T: »—1 

« T . : I" 1 
. + .. T^+e 1(1= 

A T +E T , : 11= 



Fig. 2. Approximate ?o(r) + eTi(r) temperature profiles for different [i (= —4, —1, 1) with 
5 = 2. Here, the wave number of roughness k is equal to 8, e = 0.1, (3 = 7r/4 and the inner 
boundary is s = 0.4. We can observe that the small-amplitude wavy-roughness changes the 
near-inner-wall temperature distribution significantly once /U = —4. 
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